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<~i ■ Abstract 

^ I The aim of the paper is to show that the nowadays experimen- 

. tal data from superahowed nuclear and neutron /3 decays, and from 

leptonic and semileptonic decays allow the finding of the most prob- 
able numerical form of the KM matrix, as well as the determination 
of decay constants, fp, and of various form factors f^{q^), by using 
a genuine implementation of unitarity constraints. In particular this 
approach allows the determination of semileptonic form factors that is 
illustrated on the existing data oi D ^ ttIv and D Klv decays. 



1 Introduction 

Within the Standard Model (SM) the flavour physics is encoded by the 
Kobayashi-Maskawa (KM) matrix, [1], supposed to be unitary, matrix that 
describes the quark flavor mixing through four independent parameters: 
three mixing angles, ij = 12, 13, 23, and one CP-violating phase. 
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6. By consequence the experimental determination of KM matrix entries 
is essential for the validation of the SM, and for detection of new physics 
beyond it. 

However the determination of KM entries rises a few problems. The first 
one come from theory, namely the mixing angles are not invariant quanti- 
ties, their numerical values depend on the original KM form, [T], or on the 
present day form, [2], which is not rephasing invariant, [3]. This shortcom- 
ing becomes harmless if one follows Jarlskog's approach. Starting with her 
first papers on KM matrix, she proposed the determination of the quark 
mixing matrix in terms of directly measurable quantities, and, in the same 
time, invariant quantities. In this context an invariant quantity is one whose 
numerical value does not depend on the KM matrix form, or of its rephais- 
ing invariance. Jarlskog provided two such invariants: the KM moduli, and 
the celebrated J invariant, [5j. After few years it was realized that all the 
measurable quantities of the quark mixing matrix are expressible in terms 
of four independent KM matrix moduli, [6] and [7j. In these papers it was 
also shown that areas of all unitarity triangles are equivalent, and numer- 
ically equal to half of the J invariant. As a conclusion we can state that 
the numerical values for all the measured invariant quantities should be the 
same irrespective of the physical processes where they are involved. 

The second difficulty comes from the experimental side. If one uses 
the KM moduli set as independent parameters, these ones are not directly 
measured by experimenters. In the simplest case, that of leptonic decays, 
they measure branching ratios and provide numbers for products of the form 
\Ugqi\fp, where Uqgi is the corresponding KM matrix element, and fp is the 
decay constant. For meson semileptonic decays the physical observable is the 
differential decay rate, dT/dq^, which up to known factors, is proportional 
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to \Uqqi f+{q'^)\'^ , where q denotes the transferred momentum between initial 
and final mesons. In the last case the experimental teams usually provide 
numerical values for products of the form \Uqqi f-^{0)\, where /+(0) is the 
semileptonic decay form factor at zero-momentum transfer. It is clear that 
from such measurements one cannot find two unknown parameters, say \Uqqi\ 
and fp. This can be done if and only if one can find independent constraints 
on KM matrix moduli; fortunately these ones are provided by unitarity. 
The main aim of the paper is to show how the unitarity property of the 
KM matrix can be transformed into a powerful tool for the determination of 
both matrix moduli and form factors directly from the experimental data. 

The unitarity constraints are presented in Sec. II where they are im- 
plemented in a x^ — form that depends only on KM matrix moduli. In Sec. 
Ill we present the decay formulae for superallowed 0"*" 0"*" nuclear and 
neutron P decays, and those for leptonic and semileptonic decays. They 
depend on KM matrix moduli and specific decay parameters such as decay 
constants, fp, and form factors /((?^), etc, that are implemented in an other 
X^-piece. In Sec. IV wc cite the experimental papers from which we took 
the data that we used in our fit, and in Sec. V we present the numerical 
results. The paper ends by Conclusion. 

2 Unitarity constraints 

The use of \Uij\ as independent parameters raises an important problem, 
namely the solving of the consistency problem between moduli and unitar- 
ity property, which all amounts to obtaining the necessary and sufficient 
conditions on the set of numbers \Uij\ to represent the moduli of an ex- 
act unitary matrix. After that we have to find a device for applying these 
conditions to the experimental situation where the data are known modulo 
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uncertainties. 

Both these problems have been solved, and a procedure for recovering 
KM matrix elements from error affected data was provided in [8J. These 
unitarity constraints say that the four independent parameters Sij = sin 9ij 
and cos 5 should take physical values, i.e. Sij G (0,1) and cos 5 G (—1,1), 
when they are computed via equation set: 

^ud = C^2Ci3, V^g = sf2Ci3, T^J, = S13 

Kfe = •523^13 > ^tb = ^^3023, 

^cd = •512C23 + Sl3'523Cl2 + 2S12S13S23C12C23 COS 5, 

^cs = C?2C23 + •S?2Sl3'S23 " 2si2Sl3S23Cl2C23 COS (5, (1) 

^td = •513C12C23 + Sl2'S23 - 2S12S13S23C12C23 COS 5, 

Vtl = •Si2S?3C23 + C?2'S23 + 2si2Sl3S23Cl2C23COS(5 

The above relations have been obtained by using the standard KM matrix 
form, [2j, where Vij = \Uij\, and Uij are the KM matrix entries. In the 
paper [8j it was also shown that if the independent parameters are the KM 
matrix moduli the reconstruction of a unitary matrix knowing its moduli is 
essentially unique. By consequence in the following the used independent 
parameters in our phenomenological analysis will be the Vij moduli. Al- 
though only four of them can be independent, the experimental data that 
are known only modulo uncertainties "force" us to use all the possible sets 
of four independent moduli, as it will be shown in the following. A simple 
combinatorial evaluation shows that there are 57 such sets. 

The relations ([T]) are rephaising invariant, i.e. they have the same form 
after multiplication of all the KM matrix rows and columns by arbitrary 
phases. More important is that they contain a part of the unitarity con- 
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straints. It is easily seen that all the six relations such as 



Kd + K. + Kl = i (2) 

are a consequence of the relations ([T]). The relations ([2]) are the necessary 
conditions for unitarity fulfilment, but they are not sufficient, as we show in 
the following. In fact the relations ([2]) tell us that if they are satisfied then 
there exists a physical solution for the mixing angles Sij. The relations ^ 
are not sufficient because the unitary matrices are naturally embedded in a 
larger class of matrices: the set of double stochastic matrices, [9j. 

A 3 X 3 matrix m is said to be double stochastic if its entries satisfy the 
relations 



niij 



3 3 

> 0, ^mij = l, ^mij = l (3) 
i=i j=i 



This set is a convex set, i.e. if Mi and M2 are double stochastic, then 



M = xMi + (1 - x)M2, < X < 1 



(4) 



is double stochastic. The unitary matrices are a subset of this larger set if 
we define the mij entries by the relation niij = \Uij\'^, i.e. one get relations 
similar to relation ([2]). 

For example if we choose four independent moduli, e.a: Vus = a, V^h = 
b, Vcb = c, and Vcd = d, then the following matrix 



DSi 



d^ l-c^-d^ 



b^ 

^2 



(5) 
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is a double stochastic matrix since it exactly satisfies all the six relations 

similar to relation ([2]), as it is easily checked. 

Prom any four independent moduli entering ([1]) one can get the mixing 
parameters Sij and cos 5, i.e. the four independent parameters entering the 
KM unitary matrix. For example with the above choice we find from the 
first six equations ([T]) all the three mixing parameters, Sij, and 6 as follows 

.13 = = 6, .12 = ..3 - (6) 

(1 - b^){d^l - 62) - a2) + c2(a2 + fe2(„2 ^ ^2 _ 

COS = , — , 

2a6cVl - a2 - 6Vl - P - 

and from the remaining relations three new cos 6 formulae similar to ^ . 
The above relation shows that cos 6 is an other invariant in the Jarlskog 
sense depending on four independent moduli, and the CP-violation phase 
can be measured via relations such as 

If we make use of the last four relations ([1]) we get only one solution 
for mixing parameters and cos 6. Thus depending on the chosen four inde- 
pendent moduli set the number of solutions varies between one and four. 
Because there are 57 such groups one get 165 different expressions for cos 6. 
They take the same numerical value when are computed via Eqs. ([1]), if 
and only if all the six relations similar to Eq. ^ are exactly satisfied. If the 
moduli matrix generated by four independent moduli is compatible with uni- 
tarity then cos 5 £ (—1, 1), and outside this interval when the corresponding 
matrix is not compatible. 

To see that we make use of the numerical example Vus = 2257/10^, Vub = 
359/10^ Vcd = 2256/10^ and Kfe = 415/10^. By using the DSi matrix, 
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([5]), one get the \U\ moduli as 



\U\ 



V9490466219 
105 

141 
625 

V580181 
105 



2257 
10" 

3V10526471 
W 



V10982 
2500 



359 
105 

83 
2xl0» 

V9982648619 
105 



(7) 



By using the formula from the second row of ([6]) one gets cos 6 ~ 0.64088, 
showing that the above moduli matrix, ([7]), comes from an exact unitary 
matrix. 

If we modify the previous numerical Vus value by adding to it the small 
quantity 3 x 10~^ the mixing parameters are still physical, only si2 is modi- 
fied by a very small quantity, and respectively three square root entries from 
the first two columns of ([7|), necessary for the exact fulfilment of all the six 
relations similar to In this case one gets cos 5 ~ —1.42427, which shows 
that the new moduli matrix, \U\, is not compatible with unitarity, even if 
all of its entries are positive and satisfy exactly all the six relations ([2]). If 
one computes the J invariant one finds in the above two cases 



j2 = 6.317 X 10~^°, and = -1.106 x IO- 



CS) 



Thus the physical conditions for unitarity compatibility are cos 5 G 
(— 1, 1), and > 0, respectively, and from a theoretical point of view they 
are equivalent. 

There are even bad cases; for example if instead of a and b we choose 
the parameters Vud = f, and Vts = g the corresponding double stochastic 
matrix has the form 



DS. 



■,2 



1 - C2 - c2 
■2 I f 2 



d' + P-g' 
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With the same numerical values for Vd, and Vet, and Vud = f = ^Ing^ , and Vts 



5 = 405 fij^^j 

Ki = l-c^-(i^-/^ + 5' = -^ (9) 



2 j2 f2 I 2 _ '^^'''61 



and all the other DS2 entries are positive. Thus Vub takes an imaginary 
value, and cos 6 and J take complex values. 

Similar results can be easily obtained; another example is that generated 
by the moduli, Vcd = d, Vcs = h, Vtd = Vts = 9-, whose corresponding 
double stochastic matrix has the form 



If we choose Vcd = 2252/10^ Vcs = 97345/10^ T4d = 862/10^ = 403/10^ 
one gets cos (5 ~ 0.3685 which means that the corresponding \U\ is com- 
patible with unitarity. If we subtract from Vcs value the tiny quantity 
5 X 10~® one gets cos 6 ^ —5.0234, and if we subtract 9.3 x 10~^ we get 
cos 6 —1.641 i. 

For numerical computations the use of cos 5 formulae, like ([B]), seems to 
be more efficient because of their great sensitivity to small moduli variation, 
as the above numerical computations show. 

The real physical cases are those where the central value moduli matrices, 
directly determined from data, or from a fit do not exactly satisfy all the 
six relations but only approximately; for example for a good fit the 
difference from unity could be of the order, 10"^ — 10~^, i.e. rather small 
from a phenomenological point of view. In these cases the different formulae 
for cos 6 provide different values, physical and unphysical, even if the mixing 
parameters take physical values as in the previous examples. Thus the 
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physical reality obliges us to implement the unitarity constraints 

cos5i « cos(5j, i / j, all cosJj G (—1, 1) (10) 
into a fitting device, and our choice is 

x?= E f E ^^-i) + E f E y^-^ 

j=u,c,t \i=d,s,b ) j=d,s,b \i=u,c,t j 

+ ^(cos<5(^) -cos(5(■'■))^ — 1 < cos (5^*) < 1 (11) 
that enforces all the unitarity constraints. 



3 Decay Formalism 

In this section we show how the used decay formalism for the description 
of available experimental data allows us to define a second piece of the 
X^-function by taking into account as much as possible the experimental 
information. 

For example the V^d parameter enters the description of superallowed 
0^ —7- 0^ nuclear /3 decay, and neutron /3 decay. Superallowed 0^—7-0^/3 
decay between T = 1 analog states depends uniquely on the vector part of 
the weak interaction and, according to the conserved vector current hypoth- 
esis, its experimental ft value is related to the vector coupling constant, 
which is a fundamental constant and by consequence has the same value 
for all such transitions, see |1U) . [11] . |12] and |13] . This means that the 
following relation should hold 

= 2\Gy\-\M,\- = (''^ 
where Kjijicf = 2tt^ /i ln2/(meC^)^, Gy is the vector coupling constant for 
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semi-leptonic weak interactions, and Mp is the Fermi matrix element which 
in this case is equal to \/2- The ft value that characterizes any /? transition 
depends on the total transition energy QeCi the half-life, of the parent 
state, and the branching ratio for the particular studied transition, [llj . 
However the above relation is only approximately satisfied by a restricted 
data set, and for this set one defines a "corrected" value Ft = + + 
5ns — <5c)) where 5'^ and 5ns comprise the transition-dependent part of the 
radiative correction, while 5c depends on details of the nuclear structure. In 
the above formula one takes \G'y \ = QyV^^, with gy = 1, and (|12|) is written 
in the new form 

where is the transition-independent part of the radiative corrections 
whose last estimation given in [11' is 

A]^ = (2.361 ± 0.038)% (14) 



Similarly for neutron /3 decay we make use of the formula 

V„^,a + 3A^).f?!i^ (15, 



see |14j . where Tn is the neutron mean life and \ = Qa/ dv ■ In our approach 
Vud: and A are free parameters to be found from fit. 

In SM the purely leptonic decay of a P meson, P — )• Ivi, proceeds via 
annihilation of the quark pair to a charged lepton and neutrino through 
exchange of a virtual W boson, and the branching fraction, up to radiative 
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corrections, has the form 

where Gp is the Fermi constant, Mp and rrii are the P meson and I lepton 
masses, respectively, fp is the decay constant, Vggi is the modulus of the 
corresponding KM matrix element, and rp is P lifetime. 

The next simple decays involving Vij moduli are the semileptonic de- 
cays of heavy pseudoscalar mesons, H, into lighter ones, P, whose physical 
observable is the differential decay rate, written as 

df - 192^3^3 ^ {q)\U{q)\ (17) 

where q = pn — Pp is the transferred momentum, and f+{q^) is the global 
form factor which is a combination of the two form factors generated by the 
vector part of the weak current. When the leptons are electrons, or muons 
whose masses are low compared to the mass difference, uih — mp, A(g^) is 
the usual triangle function 

A(g2) = (M^ +Ml- q^f - AMfjMl (18) 

For the decay B — )• Dlu the experimenters make use of an other variable, 
namely w = (M^ + M|, — q^) / {2MbMp,)- When the r lepton is involved 
the above formulae are a little bit more complicated, see [T5] . 

Usually, till now, the experimenters provided numerical values for prod- 
ucts of the form Vqqi /+(0), and a few for Vqqi |/+(g^)|. 

The second x^— component which takes into account the experimental 
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data has the form 




2 



(19) 



where di are the physical parameters one wants to be found from fit, di 
are the numerical values that describe the corresponding experimental data, 
while (Tj is the uncertainty associated to di. For semileptonic decays di could 
be of the form di = \f^{qf)\Vki ■ In the following our fitting x'^-function will 
be 



As it is easily seen the Vij moduli enter naturally in all the formulae 
that describe the leptonic and semileptonic decays being, in our opinion, a 
strong argument for their use as fit parameters. 

4 Experimental Data 

In our analysis we used the data on superallowed 0^ — t- 0^ nuclear /3 decays 
published in the papers [lUj-[13j. and data on the neutron lifetime from four 
papers: [16], [l7], [18], [19], for Vud determination. We also used four values 
for the /3-asymmetry parameter Aq, from papers [20], [21], [22], [23], and 
one for the electron-antineutrino correlation coefficient oq, [53], for the A 
determination. 

Vus modulus is involved in the kaon and pion leptonic and semileptonic 
decays, but also in the ratio Vus/Vudi that appears in the Marciano relation. 



2 2,2 

X =Xi+X2 



(20) 
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[25j . that we write as 



^usJ K 



udJ TT 



(1 + Cr) 



B{7r ^lytx{l))TK■mK{l 



B{K //f^(7))r^mjr(l 




(21) 



where Cr is a radiative correction stemming from both tt and K hadronic 
structures. In fact relations of the form ()2ip can be written for ah the 
leptonic and semileptonic decays, which leads to a cancellation of the exper- 
imental errors appearing on the right side. 

The last numerical values for the product f^^{0)Vus are given by KLOE 
collaboration, [26], and by FlaviaNet Working Group, [21\, and in fit we 
made use of all the (little) different /_|_(0)T4s values corresponding to the 
five channels. We also used results from [28], [29], [30], [S], [32], [33], [S], 
[35] , [36] , [37] . [38] . [39] that give only a "mean value" for the above product. 

Vub is the most poorly determined modulus although there is much ex- 
perimental information coming from decays B — t- ttIu, see [H], |42j . [43j . 
[H], [IS], [17], [IS]. In these papers the experimenters have been confronted 
with the known difficulty, getting two distinct parameters, Vub and |/+'^('?^)|, 
from their product measured from experimental data. Thus they used the 
form factor lattice computations to obtain Vub values depending on g^, see 
[41j . [42j . and [48j . For fit we found three measurements for fsVub, |43j . 
|45| . [46| . and three fenomenological determinations involving ff^{0)Vub, 



Vcd and Vcs moduli enter the leptonic decays D — )• Iv, [50], [51], [52], and 
Df — Iv, [53], [M], [55], [56], [57], respectively, as well as the semileptonic 
decays D — )• vr/z^, and D — )■ X/z^, [SS], [UD], [ST]- In the last three papers one 
find Vcqf+{0) values, for q = d, and q = s, and for the first time numerical 



m 



m 
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values for the product 



cq 



G'p 



(22) 



|61j . that ahow the form factor extraction directly from data. The above 
semileptonic decays also allow the measurement of the ratio Vcdf^^{0) /Vcsf^^{0), 



see [60], [61], [62], [63], [M], [65], and give an independent determination of 
the ratio /f ^(0)//f '^(0), which in fit was considered a new independent 
parameter. 

Finally from the semileptonic decays B — )• Dlu and B — )• D*lv, [66], [67| . 
[68], [69], [10], [71], [72], [73], [71], [75], [76], [77], one find the V^b, Q{1) and 
J-"(l) parameters. 



5 Numerical results 

Data from the above cited papers were used to define xi, the second com- 
ponent of full x^, which has a parabolic form in V^j. The first component 
Xi5 piP- that contain all unitarity constraints, has a parabolic part, and 
one that is highly nonlinear in all Vij. Thus we had to test the stability 
of the expected physical values against the strong non-linearity implied by 
unitarity. Eventually the chosen method was to modify all the measured 
central values in the same sense, plus and minus, respectively, proportional 
to their corresponding uncertainties. 

An important assumption included in our approach was: the numerical 
Vij moduli values must be the same irrespective of the physical processes 
where they are implied. Accordingly the other parameters, such as decay 
constants, /p, form factors, /+(0), A, etc., that parametrize each given ex- 
periment, have been considered as independent parameters to be obtained 
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from fit, by applying the usual technique to obtain their mean values and 
uncertainties. 

The stability tests provided sets of different moduli matrices that have 
been used to obtain a mean value matrix and its corresponding error matrix. 
The mean and uncertainty matrices have been computed by embedding the 
unitary matrices into the double stochastic matrix set, see papers [8j and 

m- 

The central values and uncertainties of data used in fit are those pub- 
lished in the above cited papers, and we combined the statistical and sys- 
tematic uncertainties in quadrature when the experimenters provided both 
of them. The numerical values obtained from fit for the decay constants, 
/tti Ik, Ib, Id, Idsj the semileptonic form factors /+(0), A^, A, as well as 
for the ratios fK/fnand, /f ^(0)//f ''(0) are given in Table 1. All of them 
are in the expected range although many have big uncertainties. 

Our approach allows a "fine structure analysis" of all experiments mea- 
suring one definite quantity, such as Apj, or f^^{0), etc. For example 
KLOE collaboration data, [26j, and FlaviaNet Working Group data, |27j . 
on f^{0)Vus lead to 

/+ "(0)xLO£; = 950.38 ± 5.56 and f^^^mFiavia = 955.06 ± 4.31 (23) 

respectively, whose central values are a little bit different, but compatible 
between them at la, and all together provide 

f+H0)KLOE+Flavia = 952.72 ± 5.30 (24) 

'^(0) from Table 1 has a precision of 1%, and the above value obtained 
from ten measurements has a precision of 0.56%. 
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Table 1. Numerical values for decay constants fp and 
form factors /+(0) in MeV units, and and A 



Parameter 


Central Value 


Uncertainty 


u 


131.131 


1.522 


fK 


154.97 


2.17 


fx/ flT 


1.1818 


0.0042 


fB 


222.8 


25.0 


fD 


207.6 


9.8 


fOs 


271.0 


18.0 


/f-(0) 


955.34 


9.27 


/f-(0) 


214.9 


13.4 




653.2 


19.1 


/f^(0) 


751.8 


10.4 


/f^'(o)//^'-(o)i 


1.171 


0.049 


^(1) 


957.5 


57.7 




1,125.3 


40.7 




2.373 


0.096 


A 


-1.2686 


0.0057 




0.002 


0.0001 



1 Data on Vcdf+'^/ycsf+^ from papers 



In this approach one could obtain information about lepton universality 
in leptonic decays. Because the corresponding /m decay constant should 
be same for both decays M — t- /xi^ and M — t- tv, a big difference between 
them could show a possible violation. As an example we chose M = Df 
meson since lattice computations provided a number for this decay constant. 
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fu^ = 241(3), see paper [78], with a very small error. Our results are 

fDt^,u = 265.0 ± 14.0, /^+^,, = 276.0 ± 20.0 (25) 

which are consistent with lepton universality. The above two numbers to- 
gether with that from Table 1 completely disagree with that provided by 
lattice computations, being far away from theoretical prediction at 8a, 13a, 
and 10a, respectively, where, a = 3, is the lattice uncertainty. The exper- 
imental spreading is, 246.0 < < 311.0, and the minimal and maximal 
values correspond to the branching ratios obtained for B{Df — t- = 
(5.15 ± 0.63 ± 0.20 ± 1.29) x 10"^^ see [SH], and to the branching ratio 
B{Df — )• Tu) = (8.0 lb 1.3 ± 0.6)%, given by Eq.(6) in paper [57|, respec- 
tively. The new result obtained in [79j does not improve the situation. 

On the other hand if one computes the difference between each one of the 
above three values and that provided by lattice computation, divided by the 
corresponding experimental error a obtained from the fit one finds the same 
value, 1.17, which shows again the consistency of experimental data. In our 
opinion the lattice number is much underestimated. In contradistinction 
their value for /x/Ztt = 1-189(7) is not far from the fit value /x/Ztt = 
1.1818(42). 

Another unexpected result concerns the constancy, usually assumed 
in all the four papers [ID]- [13], assumption that is not confirmed by our 
analysis. However our result = 2.373 ± 0.096 obtained from data [TO] 
is in good concordance with the value given by relation ()14p . see Table 2, 
but our uncertainty is 2.5a higher than the theoretical one. In this case the 

spreading that results from paper [lOj is, 2.2027 < < 2.472, which 
corresponds to 7. la where a = 0.038% is the theoretical uncertainty. The 
extremal nuclei are ^^Mg and ^^Co, respectively. The spreading obtained 
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from Savard et al data, p]3], corresponds to 10. 80", and the extremal nuclei 
are "^^Rb and ^^Cl. However the difference between the central results from 
[10| and [13j is 2cr, which suggests that there is still room for numerical 
improvements. The mean values and the corresponding uncertainties for all 
data from papers |10]-|13] are given in Table 2. 



Table 2. /x^ central values and uncertainties provided by fit when using 
data from Refs. [T0l-[T3l 



Ref. do] 


m 


m 


m 


m 


2.373(96) 


2.399(112) 


3.361(196) 


2.294(136) 


2.361(38)1 



1 The "constant" A)^ and its uncertainty, from \1X\ 

The precision of fp and /+(0) determinations in Table 1 varies from 1% 
for f^^{0) to 11% for fp- More about variability of the above parameters 
could be learnt from KM moduli matrix. Our fit result for KM central 
moduli values is 



Vr 



0.974022 0.226415 0.004251 
0.226253 0.973323 0.038108 
0.009569 0.037131 0.999265 



(26) 



Vc and its associated uncertainty matrix, fry^, have been obtained by 
using the convexity property, ([!]), and ay^ has the form 



1.1 X 10"^ 1.9 X 10"^ 2.1 X 10"^ 
3.6 X 10"^ 2.7 X 10"^ 3.1 x 10"^ 
3.5 X 10-5 2.9 X 10"^ 3.( 



1.9 X 10"^ 



(27) 



The last matrix has been obtained with the help of stabilty tests. One 
such matrix is (1281) that was obtained when all the central measured values 
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have been modified with plus one tenth from the corresponding uncertainty. 
Although such a modification is highly improbable from an experimental 
point of view, it brings to light the variation direction for all parameters 
entering the fit. 



0.974021 0.226332 0.007484 
0.226114 0.973083 0.044512 
0.012432 0.043391 0.998891 



(28) 



For example the above three matrices show that V^d is precisely deter- 
mined with four digits, while Vus^ Vcd, and Vcs only with three digits. V+ 
matrix also shows the high Vub volatility, such that future data could lead 
to higher values for it than that given by Vc matrix. In fact matrices V^. 
and Vc show that fs G (123.5,222.8). 

Although our value for Vub = (4.25ib0.02) x 10^'^ is higher than that from 
PDG fit, [2], there are experimental determinations which are compatible 
with it, one example being Ref. [H], whose value is Vub = (4.1 it 0.2^^ it 



0.2 



+0.6 
syst _04.FF- 



X 10" 



The new data from D leptonic, |50j-|52j. and semileptonic decays, [59j - 
|61j . combined with the new data from B — t- Dli' and B — t- D*!^ |66j-|70j. 
changed the KM moduli values, in particular those from the last row and 
column, see [2], p. 150. 

Our approach allows the form factor determination. The paper |61j 
provided for the first time results on Vcq\f+{q'^)\, q = d, s, and in Table 3 
are given our determinations. 

A big step forward will be the measurement of q'^ dependence for prod- 
ucts of the form \f+{q'^)Uqb\, where q = u, c, for B — )• ttIv, B — t- Dlv and 
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B D*lu decays, similar to that done for D semileptonic decays. 
Table 3. form factors from D ^ irlv and D ^ Klv decays 



Bin 


fliq^) (TrOe+z/e) 




1 


0.6895 ± 0.0576 


0.7072 ± 0.0355 


2 


0.7514 ± 0.0709 


0.7735 ± 0.0400 


3 


0.8442 ± 0.0841 


0.7956 ± 0.0488 


4 


0.8928 ± 0.0974 


0.9812 ± 0.0532 


5 


1.1005 ± 0.115 


1.017 ± 0.065 


6 


1.3083 ± 0.151 


1.101 ± 0.084 


7 


1.5779 ± 0.221 


1.635 ± 0.111 


8 




1.759 ± 0.2358 


9 




2.024 ± 0.301 








1 


0.7798 ± 0.0136 


0.7808 ± 0.0208 


2 


0.8281 ± 0.0146 


0.8096 ± 0.0239 


3 


0.8435 ± 0.0156 


0.8466 ± 0.0259 


4 


0.9113 ± 0.0198 


0.8856 ± 0.0289 


5 


0.9945 ± 0.0229 


0.9175 ± 0.0331 


6 


1.007 ± 0.027 


1.024 ± 0.039 


7 


1.128 ± 0.033 


1.149 ± 0.048 


8 


1.212 ± 0.042 


1.090 ± 0.062 


9 


1.303 ± 0.066 


1.233 ± 0.092 


10 


1.561 ± 0.165 


1.476 ± 0.210 



The simplest case is that of \f+{q^)Uui,\ because there are measured data 
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which in principle could be transformed in values for the product \f+{q'^)Uub\- 
Such measurements will allow a more precise determination of the moduli 

Vub and Vcb, and, by consequence, they will provide better values for all KM 

matrix moduli from the the first two rows, and, perhaps, the first hints for 

new physics beyond SM, if any. 

The Vc matrix, (|28|) . provides numerical values for 6 and the angles of 

the standard unitarity triangle, as follows 

6 = (89.96 ± 0.36)°, a = (64.59 ± 0.27)°, 

(30) 

7 = (89.98 ± 0.06)°, (3 = (25.49 ± 0.28)° 

The 6 value could be interpreted as a maximal violation of CP symmetry 
because sin 5 w 1. A surprising change is the shape of the standard unitarity 
triangle that is now a right triangle, since 7 w 90°. The Jarlskog invariant 
is J= (3.567 ±0.007) x IQ-^. 

Similar results are obtained from V+ matrix even if Vuu is almost twice 
bigger than that from Vc, they are 

5 = (90.0 ± 0.2)°, a = (54.1 ± 0.1)°, 

(31) 

7 = (90.0 ± 0.4)°, j3 = (36.0 ± 0.2)° 

The above results show that 5 and 7 angles are independent of Vub variation, 
while a and /3 are moderately dependent. Thus experimental results on the 
product \ f+{q'^)Uub\ could lead to a better determination of all the angles of 
the standard unitarity triangle. 

The small angles uncertainties show that Vc and V+ matrices are highly 
compatible with unitarity constraints and all the 165 cos 5 formulae provide 
very close each other values for all of them. 

A comparison with lattice computations from [50] shows that f!^'^{0) and 
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/+ (0), as well as their ratio are in good agreement with the exprimental 
values obtained by us, and our uncertainties are smaller. 

6 Conclusion 

In this paper we presented a phenomenological tool that allows determi- 
nation of the KM moduli, semileptonic form factors and decay constants 
directly from experimental data. It is based on a new implementation of 
unitarity constraints that makes use of KM matrix moduli as fit param- 
eters. These constraints are strong enough and give a consistent picture 
of nowadays flavour physics, and until now provide no signals for physics 
beyond the SM. 

A feature of our tool is that all the measurable parameters are each other 
strongly correlated, a little modification of one of them propagates to all the 
other parameters, property which is a consequence of unitarity constraints. 

The new data from the D leptonic and semileptonic decays, and those 
from B —7- D{D*)lv led to a significant change of KM moduli and of standard 
unitarity triangle shape. A crucial step forward would be the measurement 
of -B — >• Tilv form factors in bins of 0.5 GeV^, similar to that done in [6lj for 
D —7- Tilu and D — t- Klv, that will allow a better Vub determination, and a 
stabilization of moduli values entering the first two rows. 

As a final conclusion one can say that in contradistinction to experiments 
from other physics branches, those from semileptonic fiavour physics have 
to be considered as a single experiment because only using all of them, they 
could provide numbers for all the physical quantities: form factors, decay 
constants, Vij moduli, etc. 
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